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Abstract 

Inelastic rescattering effects in B decays into a pair PP of pseudoscalar 
mesons [PP = tttt or KK) are discussed. It is assumed tliat B meson de- 
cays tlirough a short-distance tree-diagram process into two jet-like states 
composed of low-mass resonances M1M2 which rescatter into PP. The 
rescattering of resonance pair M1M2 into the final PP state is assumed to 
proceed through a Regge flavour exchange. Since such processes constitute 
a fraction, diminishing with increasing energy, of total inelastic PP scat- 
tering, the inelastic rescattering contribution should die out for ^ 00. 
At TUB = 5.2 GeV, however, explicit estimates show that rescattering 
corrections could be substantial, leading to long-distance corrections to 
B^ — > KK, comparable to short-distance penguin contributions. 
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Studies of CP-violation in B decays must involve final state interaction (FSI) 
effects. Unfortunately, a reliable estimate of such effects is very hard to achieve. 
In the analyses of B ^ PP decays [P - pseudoscalar meson) only some intermedi- 
ate states, believed to provide nonnegligible contributions, are usually taken into 
account. Most often studies are restricted to the case of elastic or quasi-elastic 
rescattering. 

1. Quasi- elastic rescattering 
Consider quasi-elastic rescattering in B^ decays into tt+tt" and vr^vr^. The short- 
distance (SD) amplitudes are 

wiB'^n+n-) = -^(T + P) 

w(5°-.7rV) = (1) 

when only the dominant T (tree) and P (penguin) amplitudes are taken into 
account (Fig.l). Since isospin is a good symmetry of strong interactions, rescat- 
tering is diagonal in the basis of a definite total isospin. Consequently, one expects 
that the amplitudes W of definite final isospin I are modified by Watson phases 
6i only: 

WiB^inn)o) = --^(T + ^P) exp(z5o) 

WiB^innh) = --^Texp(z52) (2) 

The physical assumption that goes into Eq.(0) is that the probabilities of SD 
decays are not changed by such a long-distance (LD) effect as rescattering. When 
decay amplitude into a (eg.) tt+tt^ state is extracted from Eqs(H) one obtains 

W{B^ 7i+n-) = --^ exp(i52) (t + ^r(exp(A) - 1) + Pexp(iA)^ (3) 

where A = 6q — 62- The appearance of the second term on the r.h.s. of Eq.(^ 
indicates that in general final-state interactions cannot be properly described 
if one multiplies short-distance tree and penguin amplitudes in Eq.(|l]) by two 
different phase parameters [|l], |2|. The only exception is when the two phases 
6q and 62 are equal, which occurs if only elastic vrvr rescattering is allowed: the 
Pomeron does not distinguish between vr+vr" and vr^vr^, or (7r7r)o and (vr7r)2. 
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2. Short- distance penguin and rescattering contributions 

The short- distance penguin amphtude P is estimated at 20% of the tree contribu- 
tion. It is of the order of A^, where A = sin 6c is one of the four parameters in the 
Wolfenstein's parametrization of the Cabibbo-Kobayashi-Maskawa (CKM) ma- 
trix. Due to the dominance by the intermediate top quark, P has the weak phase 
(3 from the Vtd element of the CKM matrix. Although A is not small enough to 
permit real distinction between terms differing by one factor of A, we observe that 
if Fig. lb is understood as a diagram representing low-energy rescattering through 
charmed (g = c) or noncharmed {q = u) intermediate states, these processes are 
also of the order of A^. For q = c, the amplitude of this "charming penguin" has 
weak phase 0. For q = u, the diagram represents quasi-elastic (vrvr —>■ tttt) and in- 
elastic (M1M2 —* tttt) rescattering. The corresponding amplitude has weak phase 
—7 of the Vub element of the CKM matrix. With one of the goals of the present 
program of CP violation studies being the determination of CP- violating phases 
a, /3, 7, it is clearly very important to know which of the penguin-like admixtures 
to tree diagrams is in fact dominant (if any) and, consequently, which weak phase 
will be probed by experiments. 

3. Simple two-channel model 

Since a large part of the (time-reversal invariant) inelastic vrvr scattering at energy 
y/s = = 5.2 GeV (or larger) goes into multiparticle final states composed 
of noncharmed mesons, one may conjecture that inelastic rescattering processes 
shown in Fig. lb (with q = u) will be important. A very simplified two-channel 
model of what should be expected of inelastic rescattering was presented in Q . 
In this model there are two states: = Ivrvr), and I/2) representing "everything 
else" that tttt might scatter into. The most general 2x2 unitary S matrix is 



S 



cos 29 i sin 29 
i sin 29 cos 29 



(4) 



with the top left element describing tttt tttt scattering. 



If one accepts that final state interactions cannot modify the probability of the 
original SD weak decay, it follows that vector W representing the FSI-corrected 
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amplitudes is related to vector w of the original SD amplitudes through 

W = S^/^w ^ (1 + ^T) w (5) 

(where the last equality follows if S is close to 1, ie. if rescattering is assumed 
small). The appearance of the square root of the S matrix is related to Watson's 
theorem: in the basis of S-matrix eigenstates |z/), the above equation reduces to 
= e^^^Wiy, ie. the condition of unchanged probability {\W^\ = \w,y\) admits 
Watson phases only. 

With 

gl/2 

one immediately obtains 

Wi =wicose + i W2sme =\Wi\e''' (7) 
Since from nN data extrapolated to nn one gets cos2^^ ^ 0.6 — 0.7, one finds 

tan 5i = — tan ^ ^ 0.45— (8) 

W2 W2 

If B meson decays with similar strengths into the Ivrvr) and I/2) states, one cal- 
culates 61 ~ 20° — 25°. In principle, therefore, the effect of inelastic rescattering 
may be large. 

4- Inelastic rescattering at high energy 
In the previous example all states produced in the first stage of the decay could be 
rescattered into tttt. This is an overpessimistic situation. Consider a little more 
general case with three states = Ivrvr), I/2), and I/3), into which vrvr scattering 
may go. It may happen that B meson does not decay into state I/3). If vrvr 
inelastic scattering is dominated by |vrvr) I/3) transitions with |vrvr) — > I/2) 
small or negligible, then the total inelastic rescattering contribution corresponds 
to \B) — s> 1/2) I vrvr) and must be fairly weak. In such a case neglecting final 
state interaction altogether might be a good approximation: elastic FSIs result 
in the multiplication of SD amplitudes by the same phase for all isospin- (or 
SU(3)-symmetry-) related final states. 



cos 9 i sin 6 
i sin 6 cos 6 



(6) 
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In order to know whether neglecting inelastic rescattering is or is not justified, 
we have to estimate the contribution from states of type I/2) (ie. those into which 
B decays) in the unitarity relation for the / = partial wave in vrvr scattering 



Inelastic production of particles in high energy vrvr scattering is expected to 
proceed after the projectiles have exchanged a gluon. Colour separation is pre- 
vented by the production of many qq pairs whose number increases with energy. 
Colour singlet quark-antiquark pairs materialize later as resonances whose mo- 
menta are approximately parallel to the axis of tttt collision in c.m.s. Thus, a 
typical inelastic process may be represented by the diagram shown in Fig. 2a, in 
which the ordering of quark lines from top to bottom may be roughly correlated 
with the ordering in rapidity of produced resonances. Any transfer of fiavour 
quantum numbers over large distances in rapidity is suppressed. This corre- 
sponds to suppression of fiavour-exchange in quasi-two-body production processes 
TTTT M1M2 at high energies (the time- reversed process of this type is shown 
on the r.h.s. of the diagram in Fig. lb with q = u). It is through the presence 
of many free quark-antiquark pairs of unconstrained fiavour that the combinato- 
rial factors beat suppressions of fiavour transfer and, through the production of 
more and more resonances with increasing energy, they maintain approximately 
constant inelastic cross section. 

As a particular realisation of these ideas we may consider the multiRegge 
model of resonance production discussed over twenty years ago ||^. In this model 
inelastic vrvr collisions produce a variable number of resonances ordered along the 
rapidity axis with Regge exchanges in between the resonances. As energy (m^ 
in our case) increases, the average number of Regge exchanges in question and 
the average number of resonances produced must increase, so that no fiavour 
exchange over a large distance in rapidity space occurs. Thus, at high energy 
one expects the production of two-resonance states I/2) to die out. The same 
will happen to the production of three-resonance states, and so on, while the 
majority of the inelastic cross section will shift to the production of fc-resonance 
states with average k increasing with energy. This situation might be described 




(9) 



k 
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by replacing I/3) with a set of /c-resonance states \fk)- (In fact, each of states 
l/fc) should be labelled with an additional index to distinguish between states 
consisting of different k resonances.) A first thought is that states \fk) for k > 3 
cannot contribute to rescattering in B decays. However, since resonances may 
decay to lighter ones, interference terms may occur. We will argue, however, that 
such terms should be small and vanishing at high energy. 

As an example, let us consider the case when is large and the heavy 
resonance Mi from state I/2) decays into three lighter resonances, which later 
rescatter through I/4) — Itttt) (the argument is better when the number of res- 
onances into which resonance Mi decays is larger, ie. when and are 
larger). The rescattering part of the process (the diagram on the r.h.s. of Fig. 2b) 
cannot contain quark lines connecting outgoing pions: such processes are sup- 
pressed at high energies. Thus, there must be quark- ant iquark pairs produced 
out of the vacuum (eg. the line connecting the first and the third topmost res- 
onances). When one takes the square of the amplitude corresponding to the 
diagram on the r.h.s. of Fig.2b to evaluate its contribution to inelastic tttt scat- 
tering, one obtains amplitude represented by a diagram with free quark loops of 
unconstrained fiavour. These loops lead to a large combinatorial factor mentioned 
earlier. However, when the left- and right-hand sides of the diagram of Fig. 2b 
are combined to represent the decay followed by rescattering, the combinatorial 
factor is dramatically reduced: there are no quark lines of unconstrained fiavour. 
Within the language of ref.^], the resulting combined diagrams are even more 
suppressed: namely those parts of a diagram in which quark lines of an interme- 
diate resonance "cross" (the second resonance from the top in Fig. 2b) represent 
cancellations between resonances of opposite C-parities. Thus, only low-mass 
resonances Mi and M2 will contribute to the rescattering. Since with increasing 
energy inelastic vrvr scattering is dominated by channels with an ever-increasing 
number of produced resonances, it follows that inelastic rescattering in B decays 
should die out at rriB oo- It is interesting to note that the dissappearing of 
rescattering effects at ^ 00 was found also in a different framework 0. 

The mismatch between the left- and right-hand sides of the diagram in Fig. 2b 
should be attributed to the difference in the original separation of colour: receding 
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quarks of the decaying Mi resonance stretch a 3 — 3 colour string, while in the 
inelastic rescattering the dominant part comes from octet separation. The latter 
requires an exchange (diminishing with increasing energy) of a quark to reduce it 
to 3 — 3. Since vanishing of inelastic rescattering at — > C)0 is related to colour 
mismatch, this vanishing should be present in all models. 

5. Inelastic rescattering at ms = 5.2 GeV 
Although inelastic rescattering is expected to die out at high energy, one may 
ask if the value of = 5.2 GeV is sufficiently high to neglect it. In order to 
make some estimates, we turn to multiRegge models of multiparticle production 
processes In these models the distribution of the number of produced res- 
onances was evaluated as a function of energy. It turns out [0] that in the vrvr 
scattering at ^/s = 5.2 GeV, final states composed of just two resonances occur 
in 50% of cases, while three-resonance states show up in 35% cases, and so on. 
Only at much higher ^/s of the order of 20 — 30 GeV and more, the production 
of three and more resonances becomes dominant. Thus, inelastic rescattering is 
not negligible at = 5.2 GeV, and one has to estimate it explicitly. 

In order to do that, we assume that the average amplitude for the PP —>■ 
M1M2 process is approximately equal to the fiavour-exchange Regge amplitude 
in PP —>■ PP. From vrA^ data extrapolated to tttt scattering one can estimate 
that at = 5.2 GeV, the / = partial wave projection of the latter leads to 

\{PP\SReggeon\PP)\^ = 0.025 (10) 

for a typical state \PP) averaged over different SU(3)-fiavour representations of 
the PP system On the other hand, Pomeron exchange yields 

{PP\Spomeron\PP) = (tTTT | Spomeron kvr) = COS 26* ^ 0.65 (11) 

If one assumes further that all {PP\S\fk) amphtudes are approximately equal to 
{P P\S Reggeon\P P) (for average state of I/2) type this assumption is in agreement 
with phenomenological analyses [|, |T^), one obtains from Eq.(|^): 

\{PP\Spomeron\PP)\'^ + 'n'tot\{PP\SReggeon\PP)\^ = 1 (12) 

where Utot denotes the number of average channels. Using Eqs(0,|l^) one esti- 
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mates that 

ntot ~ 25. (13) 

or more (P). With 50% of the total cross section at a/s = 5.2 GeV being 

estimated as due to the production of two-resonance states, one finds that the 
number of two-resonance states is of the order of 

n2 ^ 12 (14) 



One may also try to directly estimate the number of M1M2 states produced by 
the short-distance decay mechanism. The SD decay amplitude is dominated by 
the tree-diagram contribution. The M1M2 states produced through this process 
include both states composed of light pseudoscalar mesons (vrvr in the simplest 
case) as well as various resonances. From the known 6-quark SD decay probability 

v{q^) = 2(1 - q^/mlf{l + 2q^/ml) (15) 

one can see that for itib ~ 5.2 GeV even quite massive resonances Mi (3-4 GeV) 
could be formed. The M2 spectrum extends to 1112 ~ 2.0 or 2.5 GeV . The 
average mi mass is of the order of 2.5 GeV , while m2 is a little smaller, around 
1.4 GeV . The number of M1M2 states produced in the above range of masses 
may be estimated on the basis of the ISGW2 approach [|11|. One obtains 0] 
around 10-20 different M1M2 states, in agreement with Eg . ([l^) . 



In order to perform more realistic estimates of rescattering effects and the way 
they might affect CP violation phenomena, it is important to know the phases of 
M1M2 — > PP amplitudes. These may be estimated by analysing Regge behaviour 
of such amplitudes in dual models One finds that for the FSI amplitudes of 
the r.h.s. of Fig. lb the phase factor arises from 

^] (16) 

where a{t) is the Regge trajectory. If the exchanged Reggeon just interchanges 
two quarks in the s-channel, the phase is real. After projecting Regge amplitudes 
onto the / = partial wave, the dependence on resonance masses seen in Eq.(|16D 
gives rise to a phase factor depending on the product (mim2)^. For low Mi,M2 
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masses, wherefrom the dominant part of rescattering is expected to come, this 
phase turns out to range from about —40° at (mim2)^ = 1 GeV^ to +20° at 
(mim2) ~ 7 GeV"^, where Regge approximation starts to become questionable 
(approximation of Eq.([T6|) requires values of s/(mim2)^ larger than 4 GeV~^ or 
so). 

Contributions from rescattering through different M1M2 states add up ac- 
cording to the distribution of (mim2)^ produced in SD decay. Using Eq.(^ and 
an appropriately chosen (mim2)^ distribution, the combined effect of rescattering 
through all intermediate inelastic states can be estimated from 

1 i 
{PP\W\B) = {PP\1- -Im T\PP){PP\w\B)+-J2{PP\T\f2j){f2MB) (17) 

i=i 

where the sum over j from 1 to 77-2 in the second term corresponds to the summa- 
tion over different values of (mim2)^ starting from its smallest value and according 
to its distribution. The first term in Eq. (p!7D describes SD decay amplitude with 
corrections due to elastic rescattering. Since the phases of amplitudes {PP\T\f2j) 
are not far from zero one expects that the phase of the total contribution from 
inelastic rescattering in B decays is around 90°. 

Numerical estimates of rescattering in B ^ tttt, KK decays, based on Eq . (p^Tf) , 
lead to significant corrections to SD estimates (see [Q). As an example, in Fig. 3 
the dependence on the number n of inelastic two-resonance channels included (ie. 
when n2 in Eq. (^) is replaced by variable n) is shown for the absolute size of the 
(zero isospin) W{{KK)q) amplitude. The amplitude is given in the units of the 
initial input tree amplitude T. In pure short-distance approach, this amplitude 
is equal to W{{KK)q) = Psnf^ ~ 0.1 T since the relevant penguin amplitude is 
usually estimated at Psd ~ 0.2 T. As shown in Fig. 3, the rescattering-induced 
contribution to W{{KK)q) becomes equal to the SD contribution already at 
n ^ 5. 

One has to conclude that although for ^ 00 one may expect that the 
effects of inelastic rescattering in B decays vanish, for the physical value of 
niB = 5.2 GeV these effects should be still important. Consequently, unless 
some unknown reasons justifying the neglect of rescattering exist, the extraction 
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of CP-violating parameters from nonleptonic B decays will be most probably very 
difficult. One may still hope, however, that when enough data on many different 
B decay channels are accumulated, one will be able to remove the effects of FSIs, 
thus reaching the parameters of the quark level. 
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FIGURE CAPTIONS 



Fig.l 

Dominant diagrams for B decay: (a) tree T and (b) penguin P 
Fig.2 

(a) Multiparticle production at high energy 

(b) Suppression mechanism for S ^ I/4) — > PP 

Fig.3 

Dependence of FSI-induced effects on the number of intermediate channels in- 
cluded for B — > {KK)i=o decays 
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